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Abstract 

There has recently been much interest in analytical computations of jet mass 
1 r-| ' distributions with and without vetos on additional jet activity [1-6]. An 

P-l important issue affecting such calculations, particularly at next-to-leading 

Q-i logarithmic (NLL) accuracy, is that of non-global logarithms as well as loga- 

rithms induced by jet definition, as we pointed out in an earlier work [3]. In 
this paper, we extend our previous calculations by independently deriving the 
full jet-radius analytical form of non-global logarithms, in the anti-fct jet al- 
gorithm. Employing the small-jet radius approximation, we also compute, at 
fixed-order, the effect of jet clustering on both C\ and CfCa colour channels. 
Our findings for the CpC^ channel confirm earlier analytical calculations of 
non-global logarithms in soft-collinear effective theory [5]. Moreover, all of 
our results, as well as those of [3], are compared to the output of the numerical 
program EVENT2. We find good agreement between analytical and numerical 
results both with and without final state clustering. 
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1 Introduction 

Event and jet shape variables have long served as excellent tools for testing QCD and 
improving the understanding of its properties (for a review, see [7]). Event/jet shape 
distributions have been used to extract some prominent parameters in QCD including 
the strong coupling and the quark-gluon colour ratio [8]. Due to the fact that shape 
variables are, by construction, linear in momentum, they exhibit a strong sensitivity to 
non-perturbative (NP) effects [7,9]. They have thus been exploited to gain a better 
analytical insight into this QCD domain [7,10]. Furthermore, jet shapes have been used 
not only to study the jet structure of hadronic final states, including jet multiplicities, jet 
rates and jet profiles (Ref. [11] and references therein), but also the subjet structure, or 
substructure, of the jets themselves (for a recent example, see [1]). The latter subject has 
received significant attention in recent years, particularly in the area of boosted objects 
with the aim to separate the decay products of Beyond Standard Model (BSM) particles 
from QCD background at LHC (for a review, see [12]). 
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Although shape variables are, by construction, Infrared and Collinear (IRC) safe, 
fixed-order perturbative (PT) calculations break down in regions of phase space where 
the shape variable is small. These regions correspond to gluon emissions that are soft 
and/or collinear to hard legs and lead to the appearance of large logs that spoil the PT 
expansion of the shape distribution [11] (and references therein). While measured shape 
distributions have a peak near small values of the shape variable and then go to zero, 
fixed-order analytical distributions diverge. To deal away with these divergences and 
successfully reproduce the experimentally-seen behaviour, one ought to either perform an 
all-orders resummation of the large logs, matched to fixed-order result, or rely on Monte 
Carlo event generators. We are concerned, in the present paper, with the resummation 
method as it paves the way for a better understanding of QCD dynamics including the 
process of multiple gluon radiation. The general form of resummed distributions for 
observables that have the property of exponentiation can be cast as [11] 

E(v) = C(a s ) exp [Lg 1 {a s L) + g 2 (a s L) + a s g 3 (a s L) H ] + D(v) (1.1) 

where L = ln(l/v), C(a s ) is an expansion in a s with constant coefficients that can be 
inferred from fixed-order calculations and D(v) collects terms that are proportional to 
powers of the shape variable v. The function gi resums all the leading logs (LL) a"L n+1 , 
while #2 resums the next-to-leading logs (NLL) a™L n and so on. 

There are two types of jet shape observables 2 : global and non-global [13]. Global 
observables are shape variables that are sufficiently inclusive over the whole final state 
phase space. The resummation of such variables, e.g, thrust, heavy jet mass and broad- 
ening, up to NLL accuracy have long been performed [14,15]. The resultant resummed 
distributions were then matched with NLO fixed-order results for a better agreement with 
measurements over a wide range of values of the shape variable [14, 16]. In the recent past, 
the NNLL + NLO distribution has been obtained for energy-energy correlation [17], as 
well as NNLL + NNLO [18] for the thrust distribution [19], both in e + e~ annihilation pro- 
cesses in QCD. Within the framework of Soft and Collinear Effective Theory (SCET) [20], 
the N 3 LL resummation for various event/jet variables have been performed [21] and used, 
after matching to NNLO, for a precise determination of the coupling constant a s . The 
extracted value is consistent with the world average with significant improvements in the 
scale uncertainty. 

At hadron colliders, what one often measures instead is jets, which only occupy patches 
of the phase space. The corresponding jet shape variables are thus non-inclusive, or non- 
global, and the resummation becomes highly non-trivial even at NLL level. Consider, 
for example, measuring the normalised invariant mass, p, of a subset of high-pj jets in 
multijet events. A veto is applied on final state soft activity to keep the jet multiplicity 
fixed. Jets are only defined through a jet algorithm, which generally depends on some 
parameters such as the jet size R [22]. We are thus faced with a multi-scale (p, hard 
scale, veto, jet size) problem where potentially large logs in the ratios of these scales 
appear. In addition to the Sudakov leading logs, a"ln n+1 p, coming from independent 
primary gluon emissions, there are large subleading non-global logs (NGLs) of the form 
a"ln n (a/6), where a and b are two different scales, coming from secondary 3 correlated 
gluon emissions. 

2 from the point of view of our calculations in this paper. 

3 These are emissions that are not radiated off primary hard legs. 
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We argued in [3] that in the narrow well-separated jets limit, the non-global structure 
of the p distribution, at hadronic colliders, becomes much like that of e + e~ hemisphere 
jet mass [13]. This is mainly due to the fact that non-global logs arise predominantly 
near the boundaries of individual jets. We had therefore considered e + e~ dijet events 
where only one of the jets is measured while the other is left unmeasured. We found, 
in the anti-A^ algorithm [23], NGLs in the ratio pQ/2R 2 E as well as 2E /Q where E 
and Q are the veto and hard scale respectively. These logs were completely missed out 
in [1,2]. The resummation of these NGLs to all-orders had been approximated to that of 
the hemisphere mass [13] up to terms vanishing as powers of R. Furthermore, we pointed 
out, by explicitly computing the jet mass (without jet veto) distribution under clustering, 
that different jet definitions differ at NLL due to clustering-induced large logs. Here we 
compute these logs, which we refer to as clustering logs (CLs), for the jet mass with a jet 
veto distribution. 

Within the same context of e + e~ multijet events, Kelley et al. [4] (version 1) proposed 
that if one measures the masses of the two highest-energy jets, instead of a single highest- 
energy jet as done in [3], then the resulting distribution is free from NGLs. This is clearly 
not correct since the latter shape observable, which we shall refer to, following [24], as 
threshold thrust 4 , is still non-global. To clearly see this consider, for example, the 
following gluonic configuration in e + e~ dijet events at 0(a 2 ). A gluon k\ is emitted 
by hard eikonal legs into the interjet energy region, Q. k\ then emits a softer gluon hi 
into, say the quark jet region. This configuration then contributes to the quark jet mass. 
The corresponding virtual correction, whereby gluon hi is virtual, does not, however, 
contribute to the quark jet mass. Hence, upon adding the two contributions one is left 
with a real-virtual mis-cancellation resulting in logarithmic enhancement of the jet mass 
distribution. The latter is what we refer to as NGLs. The other, antiquark, jet receives 
identical enhancement. Thus the sum of the invariant masses of the two jets does indeed 
contain NGLs contribution. The latter is actually twice that of the single jet mass found 
in [3]. 

Moreover, the authors of [4] (version 1) claimed that the anti-A^ [23] and Cambridge- 
Aachen (C-A) [25] jet algorithms only differ at NNLL for the threshold thrust 5 . From our 
calculation in [3] for the jet mass, which is not -with respect to clustering- much different 
from the threshold thrust, we know that the latter statement is incorrect. Nonetheless, 
an explicit proof will be presented below. Now, what is interesting in [4] and triggers 
the current work, is that the total differential threshold thrust distribution computed in 
the C-A algorithm and which contains neither NGLs nor CLs contributions, seemed to 
somehow agree well with next-to-leading (NLO) program EVENT2 [26]. 

In this paper we shall shed some light on the result of [4] by considering the individual 
colour, Cp, CfCa and CpT^n^ contributions to the total differential distribution as well 
as the effect of C-A clustering. We show that at 0(a 2 ) both NGLs and CLs are present 
and that the above agreement with EVENT2 is, on one side merely accidental 6 , and on 

This name is more appropriate at hadron colliders where at threshold the final state jets are back- 
to-back and there is no beam remnant [24] . 
5 This claim has been removed from version 2. 

6 As we shall see in sec. 5, while individual colour contributions do not agree with EVENT2 their sum 
does, but only in the shape variable range and for the jet-radius considered in [4]. Outside the latter 
range or for other smaller jet-radii they do not agree. 
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the other side due to the fact that the interval of the threshold thrust considered in [4] 
does not correspond to the asymptotic region where large logs are expected to dominate. 
The current work may be regarded as an extension to [3]. It includes: (a) computing 
the full R dependence of the leading NGLs coefficient in the anti-fct, (b) computing the 
small R approximation of the latter as well as the leading CLs coefficient in the C-A 
algorithm and (c) checking our findings, as well as those of [3], against EVENT2. It turns 
out, from the latter comparison, that the above approximation is actually valid for quite 
large values of R. 

While the current paper was in preparation, a paper by Hornig et al. [5] appeared in 
arXiv which studied NGLs in various jet algorithms, including anti-fct and C-A, within 
SCET. On the same day, Kelley et al. published version 2 of [4] in which they realised that 
this distribution is not actually free of NGLs and computed the corresponding coefficient 
in the anti-/c t algorithm. Our findings on NGLs, which were independently derived using 
a different approach to both papers, confirm the results of both SCET groups. Clustering 
effects on primary emission sector are unique to this paper. 

The organisation of this paper is as follows. In sec. 2 we compute the full logarithmic 
part of the LO threshold thrust distribution. We then consider, in sec. 3, the fixed-order 
NLO distribution in the eikonal limit and compute the NGLs coefficient, in both anti-fc t 
and C-A jet algorithms. In the same section we derive an expression for the CLs' first term 
as well. Note that our calculations for the C-A algorithm are performed in the small R 
limit. Sec. 4 is devoted to LL resummation of our jet shape including an exponentiation of 
the NGLs' and CLs' fixed-order terms. The latter exponentiation suffices for our purpose 
in this paper, which is to compare the analytical distribution with EVENT2 at NLO. It 
also provides a rough estimate of the size and impact of NGLs and CLs on the total 
resummed distribution. In appendix C, the corresponding resummation in SCET [4,24,27] 
is presented. Numerical distributions of the threshold thrust obtained using the program 
EVENT2 are compared against analytical results and the findings discussed in sec. 5. In 
light of this discussion, we draw our main conclusions in sec. 6. 

2 Fixed— order calculations: G(a s ) 

After briefly reviewing the definition of the threshold thrust observable, or simply the jet 
mass with a jet veto, presented in [4,24], a general formula for sequential recombination 
jet algorithms is presented. We then move on to compute the LO integrated distribution 
of this shape variable. At this order, all jet algorithms are identical. Note that partons 
(quarks and gluons) are assumed on-mass shell throughout. 

2.1 Observable and jet algorithms definitions 

Consider e + e~ annihilation into multijet events. First, cluster events into jets of size 
(radius) R with a jet algorithm. After clustering, label the momenta of the two hardest 
jets pr and pi and the energy of the third hardest jet E 3 . The threshold thrust is then 
given by the sum of the two leading jets' masses after events with E 3 > E are vetoed [4], 



pa and pl are the jet mass fractions for the two leading jets respectively. We have shown 
in [3] that the single jet mass fraction, p, is a non-global shape variable. Thus te must 
obviously be a non-global variable too. 

A general form of sequential recombination algorithms at hadron colliders is presented 
in [22]. The adopted version for e + e~ machines may be summarised as follows [22]: 
Starting with a list of final state pseudojets with momenta Pi 7 , energies E t and angles 
9i w.r.t. cm frame, define the distances 

^ = min(^ > Ef) 2(1 - j g >S ^ ) , d iB = E?, (2.2) 

where p can be any (positive or negative) continuous number. At a given stage of clus- 
tering, if the smallest distance is dij then % and j are recombined together. Otherwise 
if the smallest distance is d iB then % is declared as a jet and removed from the list of 
pseudojets. Repeat until no pseudojets are left. The recombination scheme we adopt 
here is the ^-scheme, in which pairs (ij) are recombined by adding up their 4-momenta. 
Two pseudojets, % and j, are merged together if 

2(1 -cosfly) < R 2 . (2.3) 

The anti-fct, C-A and k t algorithms correspond, respectively, to p = —l,p = and p — 1 
in eq. (2.2). We shall only consider the first two algorithms, anti— kt and C-A in this 
paper. Calculations for the inclusive kt are identical to those for the C-A algorithm as 
shown in [3]. With regard to notation, the jet-radius in [4], which we shall denote R s , is 
given in terms of R by 

R s = R 2 /A. (2.4) 

Here we work with R s instead of R. 

To verify that the definition (2.1) is just the thrust in the threshold (dijet) limit, hence 
the name, we begin with the general formula of the thrust, 

r= 1-max ^Jj^ , (2.5) 
A l^i \Pi\ 

where the sum is over all final state 3-momenta p and the maximum is over directions 
(unit vectors) n. In the threshold limit, enforced by applying a veto on soft activity, e + e~ 
annihilates into two back-to-back jets and the thrust axis, the maximum n, coincides 
with jet directions. At LO, an emission of a single gluon, k, that is collinear to, and 
hence clustered with say, pn, produces the following contribution to the thrust 

— (1 - cos0 fep J + -g-U - cos6 kpL ) + — 

where E R ^ L ) is the energy of the hard leg Pr(l) , w the gluon's energy and we have discarded 
an 0(t 2 ) term. Recalling that the first two terms in the RHS of eq. (2.6) are just the 
mass fractions p B and pl, respectively, at LO and neglecting the third term (quadratic 
in u) one concludes that 

r ~ r Eo . (2.7) 
This relation can straightforwardly be shown to hold to all-orders. 



r ~ — — (1 - cos9 kpR ) + -^-(1 - cos9 kpL ) + — (1 - cos6> fcp J(l - cos6 kpL ), (2.6) 



7 Pi may be the momenta of individual particles or each may be the total momentum of the 
particles whose paths are contained in a small cell of solid angle about the interaction point, as recorded 
in individual towers of a hadron calorimeter. 
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2.2 LO distribution 



In [3] we computed the LO distribution of the jet mass fraction, p, in the small R (R s ) 
limit using the matrix-element squared in the eikonal approximation. In this section, we 
use the full QCD matrix-element to restore the complete R s dependence of the singular 
part of the te distribution. The general expression for the integrated and normalised 
te distribution, or equivalently the te shape fraction, is given by 



E(te ,£ ) 



E 



d£. 



d 2 a 



adr' EQ dE 3 ' 



(2.8) 



'o jo 

where a is the total e + e~ — > hadrons cross-section. The perturbative expansion of the 
shape fraction £ in terms of QCD coupling a s may be cast in the form 

£ = £ (0) + S (1) + S (2) + • • • , (2.9) 

where E^ -* refers to the Born contribution and is equal to 1. The derivation of the first 
order correction, E^, to the Born approximation is presented in appendix A. The final 
result reads 
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(2.10) 



where we have used eq. (A. 2) to change the normalisation in eq. (2.8) from a to cr . The 
reason for this change is that the matrix-element we have used in EVENT2 is normalised 
to the Born cross-section 8 . The only difference between the two normalisations at 0(a s ) 
is in the one-loop constant. If we normalised to a we would have found CV(— 5/2 + 7r 2 /3) 
instead of Cp{— 1 + 71_2 /3). The function fE (R s ) is given by 



Ie (Rs 



-2 IniL In 



r ^-+2Li 2 (i? s )-2Li 2 (l-i? s ) + ^ln r ^- + 



(I) ■ w 



Notice that eqs. (2.10) and (2.11) are identical to eqs. (1) and (2) of [4] vl and the sum 
of the a s parts of eqs. (65) and (66) in [5] provided that the jet radius in the latter, which 
we refer to as R, is related to R s by: tan 2 (.R/2) = R s /(1 — R s )- It is worthwhile to note 
that in the limit R s — > 1/2 the te distribution (2.10) reduces to the well known thrust 
distribution [28] with upper limit r < 1/3. For R s < 1/2 the threshold thrust distribution 
includes, in addition to thrust distribution, the interjet energy flow distribution [29] too, 



2vr 



-41n- 



R.t 



-In 



1 — R s V Q 



o 



Eo 
Q 



(2.12) 



Here the interjet region (rapidity gap), referred to in literature as A?7, is defined by the 
edges of the jets. Specifically, it is related to the jet-radius R s by 



Arj 



-In 



R s 



R f 



(2.13) 



8 Note that there are three sets of matrix-elements included in the program, of which only one is not 
normalised to the Born cross-section. 
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The important features of the te distribution that are of concern to the present paper 
are actually contained in the second order correction term E^ 2 \ which we address in the 
next section. 



3 Fixed— order calculations: 0(p? s ) 

We begin this section by recalling the formula of the matrix-element squared for the 
e + e~ annihilation into two gluons, e + e~ — > liPa) + q(Pb) + 9i(ki) + #2(^2) in the eikonal 
approximation. Let us first define the final state partons' 4-momenta as 



Pa 
Pb 



f(i,o,o,i) 



Q , 



1,0,0,-1), 



u>i(l, sin 61 cos 0i, sin 61 sin 0!, cos 9i), 
002(1, sin 62 cos 02, sin 6 2 sin 2 , cos 9 2 ). 



(3.1) 



where the angles 8i are w.r.t. p a direction (which lies along the z-axis) and we assume the 
energies to be strongly ordered: Q uo\ 002- This is so that one can straightforwardly 
extract the leading NGLs. Contributions from gluons with energies of the same order, 
Q 3> uji ~ Co>2; are subleading and hence beyond our control. The recoil effects are 
negligible in the former regime and are thus ignored throughout. The eikonal amplitude 
reads [11], 

S ab (h, k 2 ) = C 2 F W P + C F C A W S , (3.2) 

where Wp and Ws stand for primary and secondary emission amplitudes respectively. If 
we define the antenna function Wij(k) = 2(ij)/(ik)(kj) then the latter amplitudes are 
given by 

16 

W P = w ab (ki)wab(k 2 ) = 2 2 . 2 . 2 , (3.3) 



and 



W s 



Wabih 



oj\oj\ sin 9f sin Q\ 



[w a i(k 2 ) + w b i(k 2 ) - w ab (k 2 )} 



u(u2 sin 9( sin 9 



2 

2 L 



1 — COS 6\ COS 62 



cos 9 



12 



(3.4) 



For completeness, the two-parton phase space is given by 



&<$> 2 {kx,k 2 ) 



i=l 



sin OidOi 
2^T 



2^ 



(3.5) 



It is worth noting that the primary emission, Wp, contribution to the te distribution is 
only fully accounted for by the single-gluon exponentiation in the anti-fct algorithm case. 
If the final state is clustered with a jet algorithm other than the latter, W p integration over 
the modified phase space, due to clustering, leads to (see below) new logarithmic terms 
that escape the naive single-gluon exponentiation. On the other hand, the secondary 
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amplitude Ws contribution is completely missing from the latter Sudakov exponentiation 
in both algorithms. 

First we outline the full a 2 s structure of the te distribution up to NLL level in the 
anti-fct including the computation of the NGLs coefficient. After that, we investigate the 
effects of final state partons' clustering on both primary and secondary emissions. The 
C-A algorithm is taken study to illustrate the main points. Calculations where 

the final state is clustered with other jet algorithms should proceed in an analogous way 
to the C-A case. 

3.1 te distribution in the anti— k t algorithm 

The anti-/c( jet algorithm works, in the soft limit, like a perfect cone. That is, a soft gluon 
hi is clustered to a hard parton pj if it is within an angular distance 2\/RZ, (= R), from the 
axis defined by the momentum of the latter. This feature of the algorithm greatly sim- 
plifies both fixed-order and resummation calculations. Considering all possible angular 
distances between (k±, k 2 ) and (p a ,Pb) we compute below the corresponding contributions 
to primary and secondary pieces of the r Eo distribution. Note that we use LL and NLL 
to refer to leading and next-to-leading logs of r Eo (and not 2E /Q) in the exponent of 
the resummed distribution (discussed in sec. 4). 

3.1.1 C% term 

The LL contribution to the te distribution comes from diagrams corresponding to two- 
jet final states. That is diagrams where both real gluons, k\ and k 2 , are clustered with 
the hard partons p a and pb- Diagrams where one of the two gluons is in the interjet 
region, and hence not clustered with either hard parton, contribute at NLL level. Other 
gluonic configurations lead to contributions that are beyond our NLL control and thus not 
considered. The Cp part of the 0(ot 2 s ) te distribution may be found by expanding the 
exponential of the LO result (2.10). The full expression including the running coupling 
at two-loop in the MS will be presented in sec. 4. For the sake of comparison to the 
clustering case, we only report here the the LL term, which reads 



Next we consider the derivation of the CfCa contribution to the te distribution 
including the full jet-radius dependence. 

3.1.2 C F C A term and NGLs 

In the anti-fct algorithm the non-global logarithmic contribution to the te distribution 
is simply the sum of that of the single jet mass fraction, p, with a jet veto distribution 
studied in [3] 9 . This is in line with the near-edge nature of non-global enhancements. In 
two-jet events, the well separated 10 jets receive the latter enhancements independently 
of each other. Possible final state gluonic arrangements relevant to NGLs at second order 

9 Here we go beyond the small R s approximation assumed in [3]. 

10 such that the jet-radius is much smaller than the jets' separation; R s -C (1 — cos 9ij), where 9ij is 
the angle between jets i and j. 




(3.6) 



S 




Figure 1: Schematic representation of gluonic arrangement giving rise to NGLs. We have 
only shown the NGLs contributions to the p#-jet. Identical contributions apply to the 
PL-jet. 



are depicted in fig. 1. The all-orders resummed NGLs distribution may be written in the 
form [13] 

S(t) = l + S 2 t 2 + ■■■ = ! + J2 S nt n , (3-7) 



n=2 



with t being the evolution parameter defined in terms of the coupling a s by 

t - 
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n dk t 


2^ 


Jfcmin kt 


a s 


/ Lmax N 

In ( 7 . 


2^ 





HJ, 



(3.8) 



where the exact form of the upper and lower limits, kf 13 * and fcf m , depend on the gluonic 
configuration and the second line in (3.8) assumes a fixed coupling. To make contact 
with interjet energy flow calculations [30,31], we work in this particular section with 
hadronic variables (k t ,r],<f)) instead of e + e~ variables (E,9,(f>). The pseudo-rapidity 77 
and transverse momentum k t (both measured w.r.t. incoming beam direction) are related, 
respectively, to the angle and energy by 11 



-In ( tan - ) , E = k t cosh(?7) 



(3.9) 



Using the secondary emissions eikonal amplitude (3.4) in terms of the new variables, the 
NGLs coefficient 5*2 reads 



S, = —AC wC a / d$ (2) 



cosh^x - rj 2 ) 



cosh^! — rj 2 ) — cos( 



- 1 



(3.10) 



where the phase space measure, d& 2 \ is of the general form given in eq. (3.5) with the 
k t integrals included in the definition of t (3.8) and new restrictions coming from the jet 
shape definition. For configuration (a) in fig. 1, it reads 



d$i 2 ) 



At) 



drj] 



+00 



d V2 ^ 6 ( ln-^ 2 - - V2 )Q(E - k tl cosh(r /l )) , (3.11) 



11 Otherwise, one can redefine the partons' 4-momenta in terms of 77 and kt and use the antenna 
function expressions of Wp and Ws to rewrite them in terms of the hadronic variables. 
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where the interjet (gap) region, Arj is given in eq. (2.13). Due to boost invariance of 
rapidity variables the latter region has been centred at rj = 0. Moreover, the factor 
2 in (3.11) accounts for the p_L~jet contribution. Since neither the integrand nor the 
integral measure in eq. (3.10) depends explicitly on the azimuthal angles (fa, fa), we use 
our freedom to set fa = 0, average over fa and then perform the rapidity integration. 
The resultant expression for S 2 m configuration (a) at the limit te — > reads, 



5*2,a — —4:C f Ca 



^ + Ar/ 2 - A V In (e 2A " - l) - ±Li a (e" 2 ^) - l -U 2 (l - e 2A ") 



(3.12) 

An identical expression was found for the NGLs' coefficient in the interjet energy flow 
distribution [30] 12 . The fact that S^a is the same for te and interjet energy flow 
distributions means that the NGLs' coefficient only depends on the geometry of the 
phase space and not on the observable itself. This is of course only true in the limit 
where the jet shape variable goes to zero. The difference between the jet shape variables 
amounts only to a difference in the logarithm's argument. 

It should be understood that there are 0-function constraints on k tl and k t2 resulting 
from rapidity integrations not explicitly shown in eq. (3.12). Performing the remaining 
trivial k t integrals yields 

where a factor of 1/2 has been absorbed in S 2 ^ a (3.12). 

Now consider configuration (b) in fig. 1. Adding up the corresponding virtual correc- 
tion, one obtains the following phase space constraint 



n'''-n^jn fe -«j' (3i4) 

The phase space measure d$[ 2 ^ is analogous to d<3>i 2 ^ in (3.11) with 1 <H> 2 and the 
two O-functions in (3.11) replaced by those in eq. (3.14). The limits on rji are then 
+oo > rji > max[Ar)/2,\n(k t i/QTE )]. If we impose the constraint given in eq. (3.13), i.e, 
2E /Q 3> te /R s , then the lower limit becomes rji > \n{k t \/QrE )- The NGLs coefficient 
5*2,6 thus reads 

S 2 , b = -±C F C A / d Vl / d V2 [cothfa - Tfc) - 1] , (3.15) 

where we have averaged the eikonal amplitude W s over <p 2 and moved k t s' 0-functions 
onto the integral of the evolution parameter t^, which is given at a 2 s by 

The 5*2,6 1 2 contribution is then beyond our NLL accuracy. In fact, 52,6 vanishes in the 
limit te — > as can be seen from eq. (3.15). 



2 Our jet-radius, R s , is given in terms the parameter c, used in [30], by the relation: 1 — c = 2R S 
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The last contribution to NGLs at 0(a 2 ) comes from configuration (c) in fig. 1. Upon 
the addition of the virtual correction, one is left with the constraint 

6 (Qt Eo - k tl e~^) (k t2 e + ^ - Qt Eo ) . (3.17) 

The corresponding NGLs coefficient and evolution parameter read 



+00 



1t Eq 



2 



A?7 
2 



S 2 ,c = -4C F C A dru d V2 [cothfa - Tfy) - 1] , (3.18) 



l Qte 

(g)V(r £o e^/ 2 ), (3.19) 

Since we have assumed strong ordering, k t \ 3> kt2, then the lower limit of rji in (3.18) is 
ln(k t i/QTE )- Consequently the coefficient S 2 , c vanishes in the limit r Eo — > 0. For this 
reason, this configuration will not be considered. 

We conclude that in the regime 2E /Q ^> t Eo / R s , the only non-vanishing contribution 
to the NGLs comes from the phase space configuration (a). Other configurations, (b) and 
(c), vanish in the limit t Eq — > 0. Hence 

S 2 = S 2 , a , t = t a . (3.20) 

In fig. 4 we plot S 2 as a function of the jet-radius R s . At the asymptotic limit Arj — > +00 
(or equivalently R s — > 0) S 2 saturates at — C e Ca 2tt 2 /3. This value (or rather half of it) 
is used as an approximation to S 2 in [3]. From eq. (3.12), we can see that the correction 
to such an approximation is less than 10% for jet-radii smaller than R s ~ 0.28, which is 
equivalent to R ~ 1. Furthermore, Eq. (3.12) confirms the claim made in the same paper 
that NGLs do not get eliminated when the jet-radius approaches zero. One may naively 
expects that when the jet size shrinks down to (R s — > 0) there is no room for gluon k 2 
to be emitted into. This means that t Eo becomes inclusive and hence S 2 vanishes. To 
the contrary, S2 reaches its maximum in this limit. 
Few important points to note: 

• If we choose to order the energy scales in the ©-functions of (3.13) and (3.16) 
the opposite way, i.e, 2E /Q t Eq /R s then configuration (6) becomes leading, 
in NGLs, while the contribution from configuration (a) vanishes. That is t\ in 
eq. (3.16) becomes 

a _ ( a s\ 2 1 2 ( ^EqR s \ f t Eo 2E 



*i = ^ n^J e lf-^J' (321) 

and S 2j b = S 2 ^ a in eq. (3.12). We do not consider this regime here though. 

If, on the other hand, we do not restrict ourselves to any particular ordering of 
the scales, as it is done in Refs. [5] and [4], then both configurations (a) and (6) 
would contribute to the leading NGLs. Adding up t 2 a , in (3.13), and t 2 , in (3.21), 
the ©-functions sum up to unity and one recovers the result reported in the above 
mentioned references. Notice that it is a straightforward exercise to show that 
eq. (3.12) is equal to /ol+/or given in eq. (28) of [5] in the case where R E = Rr = R 
(= R given in sec. 2). Moreover, the coefficient /^cl gi ven i n ec L- (B2) of [4] v2 is 
related to S 2 , a by f N £ L = -8 x S 2 , a . 
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• Setting the cut-off scale E ~ te Q in t a , eq. (3.13), and tf,, eq. (3.21), would 
diminish NGLs coming from both configurations (a) and (b) and the threshold 
thrust becomes essentially a global observable. This is unlike the observation made 
in the study of the single jet with a jet veto distribution [3] where the above choice 
of E kills the NGLs near the measured jet but introduces other equally significant 
NGLs near the unmeasured jet. 

In the next subsection we recompute both C\ and CfCa contributions to the te 
distribution under the C-A clustering condition. For the CfCa term, we only focus on 
configuration (a) and do not attempt to address the subleading contributions coming 
from configurations (b) and (c). 

3.2 te distribution in the C— A algorithm 

The definition of the C-A algorithm is given in eq. (2.2) with p = 0. Unlike the anti-/c t 
algorithm, which successively merges soft gluons with the nearest hard parton, the C-A 
algorithm proceeds by successively clustering soft gluons amongst themselves. Conse- 
quently, a soft parton may in many occasions be dragged into (away from) a jet region 
and hence contributing (not contributing) to the invariant mass of the latter. The jet 
mass, and hence te , distribution is then modified. It is these modifications, due to 
soft-gluons self-clustering, that we shall address below. 

Any clustering-induced contribution to the te distribution will only arise from phase 
space configurations where the two soft gluons, ki and k 2 , are initially (that is, before 
applying the clustering) in different regions of phase space. Configurations where both 
gluons are within the same jet region, gluon k\ is in one of the two jet regions and gluon 
k 2 is in the other or both gluons are within the interjet region are not altered by clustering 
and calculations of the corresponding contributions will yield identical results to the anti- 
kt algorithm. We can therefore write the te distribution in the C-A algorithm, at 0(a^), 
as 

4 2 l A (r £o , E ) = ^ ti _ kt (r Eo ,E ) + 5^ 2 \r Eo , Eq). (3.22) 

It is the last term in eq. (3.22) that we compute in the present subsection. Starting at 
configurations with two gluons in two different regions, the jet algorithm either: 

(A) recombines the two soft gluons into a single parent gluon if the clustering condi- 
tion (2.3) is satisfied. The latter parent gluon will either be in one of the two jet 
regions or out of both of them (and hence in the interjet region). 

(B) or leaves the two gluons unclustered, if the clustering condition is not satisfied. This 
case is then identical to the anti-fct one but with a more restricted phase space. This 
restriction comes from the fact that for the two gluons to survive the clustering they 
need to be sufficiently far apart. Quantitatively, their angular separation should 
satisfy the relation 

(1 -cos #i 2 ) > 2R S . (3.23) 

Below, we examine the contributions from configurations (A) and (B) to the C\ and 
CfCa colour pieces of the te distributions. All calculations are performed in the small 
R s approximation using the e + e~ variables (u,9,(f)). 
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Figure 2: A schematic representation of a three-jet final state after applying the C-A 
algorithm on real emission along with virtual correction diagrams. The two gluons are 
clustered in the E-Scheme (see sec. 2). Identical diagrams hold for the p^-jet. 



3.2.1 C% term 

Consider the gluonic configuration in (A) where the harder gluon k\ is in the interjet 
region and the softer gluon k 2 is in the p_R-jet region. We account for the p^-jet region 
through multiplying the final result by a factor of two. Applying the C-A algorithm (2.2), 
the smallest distance is d min = d\ 2 . Hence gluon k\ pulls gluon k 2 out of the PR-jet region 
and form a third jet, as depicted in fig. 2. The latter is then vetoed to have energy less 
than Eq. The corresponding clustering angular function, in the small angles limit, reads 

e c - A (i,2) = e(0?-4i* a )e(4i*,-0f)e(0f-0? 2 ), 



q(46 2 2 cos 2 <p 2 - ei)e(ef - ar s )6(ar s - e 2 )Q e 2 ±— e cos r 

COS z (Do I V 2 



(3.24) 



Adding up the corresponding virtual corrections, where one or both of the gluons are 
virtual, one obtains the following constraint on the phase space 

B(E -ui- cu 2 ) - B(E - m) + {^e\ - r E ^ . (3.25) 

Since we are working in the strong energy-ordered regime, u\ ^> u 2 , only the last 
function survives. The new contribution to the Cp piece of the Tp distribution is then 
given by 

P2 o f Q/2 du 2 f Q/2 Au l p d0 2 f 2 ^ cos ^d6 1 f R ^ cos ^d6 2 

Go L 



J?IZl oj 2 J W2 Ul J* 2vr J 2VWa Q x J 2VWs 6 2 ' 
= 0.73C*(£)Vg). ,3.26) 

This result is identical to 13 that found in [3] for a single jet mass (without a jet veto) 
distribution. The reason for this is that the clustering requirement only affects the dis- 
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virtual corrections 



k 2 




Figure 3: A schematic representation of a two-jet final state after applying the C-A 
algorithm on real emission along with virtual correction diagrams. The two gluons are 
clustered in the E-Scheme (see sec. 2). Identical diagrams hold for the left (pl~) jet. 

tribution to which the softest gluon contributes. Which is in both cases the jet mass 
distribution. 

The second possible configuration that corresponds to case (A) is where gluon k\ is in, 
say, the p_R-jet region and the softer gluon k 2 is in the interjet region. If the two gluons 
are clustered, i.e, gluon k\ pulls in gluon k 2 , then upon adding real emission and virtual 
correction diagrams, depicted in fig. 3, one obtains the following phase space constraint 

_e { te ° - H^) 9 {% el - te ) + e{uj2 - Eo) ' (3 - 27) 

where we have assumed small angles limit and employed the LL accurate approximation 

Given the fact that u>i ^> u 2 and 9\ and 9 2 must be close to each other to be clustered, 
i.e, they should satisfy condition (3.23), then the first two 0-functions in eq. (3.27) are 
substantially suppressed and one is only left with the veto on uj 2 . Applying the C-A 
algorithm one obtains an identical clustering function to eq. (3.24). Hence the CLs' 
coefficient for this configuration is equal to C 2 given in eq. (3.26). That is C 2 = 0.73 C\. 
The evolution parameter does however change. It is now given, at 0(a 2 s ), by 




(3.29) 



This contribution is then beyond our NLL control. Note that the CLs contribution in 
eq. (3.29) is equal to what one would find for interjet energy flow distribution provided 
that the rapidity gap is defined through eq. (2.13). 

Let us now turn to case (B) where the two gluons are not merged together. If gluon k\ 
is in the interjet region and gluon k 2 is in one of the two jet regions then the corresponding 
phase space constraint reads 

e - e( Wl - Eq) [i - e c _ A (i, 2)] . (3.30) 

13 It is actually twice 
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The limits on #2~integral are then given by: min(4i? s , 2Qt Eo /uj2) > O 2 > 0. Imposing the 
constraint 2E /Q ^> te /R s , it is straightforward to see that the above constraint yields 
NNLL contribution and thus beyond our control. Similarly, the configuration where gluon 
ki is in the jet region and gluon k 2 is in the interjet region yields subleading logs. 

Hence the Cp piece of the clustering-induced correction term 8Y>( 2 \ in eq. (3.22), up 
to NLL, reads 

5^ 2 \r Eo ,E ) = C^t 2 p , (3.31) 
Next we compute the C F C A piece of SYP^ . 



3.2.2 C F C A term 

Consider the gluonic configuration (a) depicted in fig. 1. Applying the C-A clustering 
algorithm on the latter yields two possibilities. Namely the two gluons are either clustered 
or not. The former case completely cancels against virtual corrections and thus does not 
contribute to NGLs. It is when the two gluons survive the clustering, the latter case, that 
a real-virtual mismatch takes place and NGLs are induced. The corresponding evolution 
parameter is equal to t of the anti-kt case, eq. (3.20). The clustering condition is simply 
one minus that in eq. (3.24). The NGLs' coefficient can then be written, using the eikonal 
amplitude (3.4), as 

St A = S 2 + 5V { cic A (3-32) 
where 52 is given in eq. (3.20) and 

»20 2 cos<fe ^ ,2^7 ^ /-7T/3 d - 



SZ® a =8C F C A [ ™±-f 

v cos 4>2 



tt/3 



2tt 



— cos 6i cos 6 2 
1 — cos 9 
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X 



x 9 ( Rs ^-Y (3.33) 

\T Eo COS(p2 2w 2 / 

We can perform the ^i-integral analytically and then resort to numerical methods to 
evaluate the remaining 62 and 02 integrals. The result, in terms of the jet-radius R s , is 
depicted in Fig. 4. -S 1 ^^ saturates at around 0.44 x 2n 2 /3 C F C A ~ 2.92 C F C A , i.e, a 
reduction of about 55% in S 2 - This is due to the fact that for the two gluons to survive 
clustering they need to be sufficiently far apart (6* 12 > R = 2y/R~ s ). The dominant 
contribution to S 2 comes, however, from the region of phase space where the gluons are 
sufficiently close. This corresponds to the collinear region of the matrix-element; Q\ ~ 82- 
Hence the further apart the two gluons get from each other, the less (collinear) singular 
the matrix becomes and thus the smaller the value of NGLs coefficient. 

Note that the C-A coefficient S%~ A = 2 x /qr A ' wnere /or A ^ s gi ven by eq. (38) in [5], 
at least in the small jet-radius region. Noticeably, the two results coincide at both limits 
R s — > and R s — > 1/2 (equivalently R — > and R — > y/2 in [5]). In fact, the coefficient 
5 2 C " A is valid, as we shall see in sec. 5, for quite large jet-radii; up to R s ~ 0.3 (equivalent 
to R ~ 1 in [5]). 

The fixed-order NLL logarithmic structure of the t Eq distribution should by now be 
clear for both jet algorithms. In order to assess the phenomenological impact of NGLs and 
clustering requirement on the final cross-section, it is necessary to perform an all-orders 
treatment, which we do below. 
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Figure 4: Non-global coefficient 5*2 in the anti-A^ and C-A algorithms. 

4 Resummation of t# distribution 

Resummation, which is essentially the organisation of large logs arising from soft and/or 
collinear radiation to all-orders, is based on the factorisation property of the pQCD 
matrix-element squared for multiple gluon radiation. This is only true for independent 
primary emissions though. Including secondary correlated emissions, the picture dramat- 
ically changes and the resummation can only be performed at some limits, eg. large-iV c 
limit [32]. In the standard method [14,33,34], resummation is carried out in Mellin 
(Laplace) space instead of momentum space. Only at the end does one transform the 
result back to the momentum space through (inverse Mellin transform), 

£p(t £o ,£ ) = J ^- e^o J e » E ° V P { V ~\ (4.1) 

where P stands for primary emission. With regard to non-global observables, the impor- 
tant point to notice is that the resummation of NGLs is included as a factor multiplying 
the single-gluon Sudakov form factor, Sp, [13] 

£(rp ,£ ) = £p(rp ,£ ) S (t) , (4.2) 

In this section, we first consider resummation of rp distribution in events where the final 
state jets are defined in the anti-fct algorithm and, second, discuss the potential changes 
to the resummed result when the jets are defined in the C-A algorithm instead. 

4.1 Resummation with anti— kt algorithm 

As stated in the introduction and proved in sec. 2, the Tp observable is simply the sum 
of the invariant masses of the two highest-energy (or highest-p t for hadron colliders) jets. 
Therefore the Tp resummed Sudakov form factor is just double that computed in [3], for 
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a single jet mass. That is, up to NLL level we have 



exp 


-2 


[Kt Eo (t Eo )+1eK Eo (te )) 




r 


[i + 2n> TEo (T Eo )) 





exp [-1Ze (E )] . 



(4.3) 



The full derivation of (4.3) as well as the resultant expressions of the various radiators 
are presented in the small jet-radius limit in Ref. [3]. To restore the full R s dependence 
we make the replacement R 2 /p R s /{r Eo {l — R s )) such that when expanded eq. (4.3) 
reproduces at 0(a s ) the LO distribution (2.10). 

To account for the NGLs at all-orders in a s , it is necessary to consider an arbitrary 
ensemble of energy-ordered, soft wide-angle gluons that coherently radiate a softest gluon 
into the vetoed region of phase space [13] 14 . The analytical resummation of NGLs is then 
plague with mathematical problems coming from geometric and colour structure of the 
gluon ensemble. Two methods have been developed to address this issue: A numerical 
Monte Carlo evaluation [13, 30] and a non-linear evolution equation that resums single 
logs (SL) at all-orders [35]. Both methods are only valid in the large-iV c limit. In the 
latter limit and for small values of the jet-radius R s , we argued in [3] that the form of 
S(t) should be identical to that found in the hemisphere jet mass case [13]. Since in the 
present paper we are not confined to the small R s limit, we need to modify and re-run the 
Monte Carlo algorithm, presented in [13], for medium and large values of the jet-radius 
should we seek to resum the te NGLs distribution. The latter task is, however, beyond 
the scope of this paper. Here, we are only aiming at comparing the analytical results 
with fixed-order NLO program EVENT2. It suffices in this case to simply exponentiate 
the first NGLs term in eq. (3.20), 

S(t) = exp (S 2 t 2 ) . (4.4) 
The distribution (4.2) is of the generic form given in eq. (1.1). Explicitly, it reads 



E(t Eo ,E ) 




exp 



OO 71+1 



n=l m=0 



+ Aln(T£ ), (4.5) 



where Ck is the k th loop-constant, L = ln(l/rE ) and = D, which vanishes in the 
limit te —> 0. In order to determine the coefficients G nm at NLO and up to NLL,we 
need to expand the radiators, as well as the T function, in eq. (4.3) up to second order 
in the fixed coupling a s = a s (Q). The results are presented in appendix. B. Although 
we have provided the NNLL coefficient, G21 in eq. (B.l), we do not claim that it is under 
control. Nonetheless, it does capture all -R s -dependent terms 15 . The missing terms from 
G21 include: a) coefficients of L which are independent of \n(R s /(l — R s )) for all colour 
channels. These can be borrowed from thrust distribution [19,36,37]. b) Although G21 
has a sub leading NGLs term in the C e Ca colour channel, which comes solely from the 
expansion of t (eq. (3.20)), the full expression in this channel as well as in the CfTrUj 



14 In our case the vetoed region is the jet region. Due to symmetry, we can choose one jet region and 
multiply the final answer by a factor of two. 

15 as can be seen from comparison to the SCET result (C), which only contains the primary emission 
piece and is valid to NNLL. 
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channel is still missing. To properly compute the latter, one has to extend both the 
matrix-element (3.2) and the phase space to include hard emission. Such a task will 
be considered elsewhere. It is worthwhile to mention that full subleading NGLs have 
recently been computed analytically within SCET framework for the hemisphere mass 
variable [36,37] 16 . The two-loop constant C2 has also been computed for the latter 
variable as well as the thrust [19,36]. 

To make contact with SCET calculations, we provide in appendix C the full formula 
of the Sudakov form factor for the te primary distribution including determination of 
the G nm coefficients in SCET. 

Next we comment on the form of resummation when final state jets are defined in the 
C-A algorithm. 

4.2 Resummation with C— A algorithm 

With regard to primary emission piece, resumming logs induced by clustering is a cum- 
bersome but doable task. It has been performed, for example, in [38] for interjet energy 
flow distribution where final state jets are defined in the inclusive k t algorithm. The final 
result of the resummed radiator was written as an expansion in the jet-radius and the 
first four terms were determined. For secondary emissions, the resummation of NGLs has 
only been possible numerically and in the large- N c limit. It has again been carried out 
for the above mentioned energy flow distribution in [31]. We expect that analogous, to 
the interjet energy flow, analytical treatment and numerical evaluation can be achieved 
for the resummation of CLs and NGLs, respectively, for the r^ variable. We postpone 
this work to future publications. 

For the sake of comparing to EVENT2, it is sufficient to simply exponentiate the fixed- 
order terms S^ _A and C<f , just as we did with the anti-A^ algorithm case. Due to the fact 
that logarithmic contributions induced by clustering arise mainly from soft wide-angle 
gluons, we expect them -clustering-induced logs- to factorise from the primary form factor 
at all-orders. Therefore, the resummed distribution, whereby clustering is imposed on 
the final state, may be written as 



where S c A is of the form (4.4) with S2 replaced by S% A and, in analogy with the NGLs 
factor, the CLs factor reads 



In fig. 5 we plot the resummed differential distributions, dE(re , E )/(1te = (l/<7 )da/dr El 
computed from eq. (4.3) for the anti-A^ algorithm and from eq (4.6) for the C-A algo- 
rithm at different values of R s . The dependence on Eq has been discussed in [3] where the 
all-orders NGLs resummed expression was employed. There are several points to note. 
Firstly the effect of NGLs is a suppression of the total cross-section relative to the pri- 
mary result. This suppression is diminished by decreasing the value of R s . For example, 

16 Recall that for the leading NGLs the corresponding coefficient for the hemisphere mass distribution 
corresponds to setting R s = in S2.11 (3.12). 



E(r £o , Eq) = Sp(r Eo , Eq) S c ~ a (t) C p (t p ) , 



(4.6) 




(4.7) 
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at E = 0.1Q the Sudakov peak is reduced due to NGLs by about 4.02%, 3.42%, 2.62% 
and 1.35% for R s = 0.30,0.12,0.04 and 0.0025 (equivalent to R = 1.1,0.7,0.4 and 0.1) 
respectively. These values are only meant to give an idea of the effect of varying the 
jet-radius parameter on both NGLs and CLs corrections to the total cross-section, since 
we are only working with an approximation of the latter and not the full all-orders result. 
It has been shown in [31], for the interjet energy distribution, that the NGLs resummed 
factor S(t) at all-orders is much smaller (thus larger suppression of primary-only result) 
and of different shape, as a function of t, to the fixed-order exponentiated result. 

Secondly the effect of clustering is reducing the phenomenological significance of NGLs. 
This reduction becomes larger, hence the NGLs suppression on the Sudakov peak becomes 
smaller, as one moves towards smaller values of R s . For the same jet veto Eq = 0.1Q, the 
Sudakov peak is reduced by 0.62% 17 , 0.80%, 0.63% and 0.22% for R s = 0.30,0.12,0.04 
and 0.0025 respectively (values are only an estimate of the impact of clustering). Com- 
paring to the anti-fct case, we see that the effect of NGLs has been reduced by more than 
70% for R s = 0.12,0.04 and R s = 0.0025. This observation suggests that instead of re- 
summing NGLs, which is a daunting task even numerically, one should, perhaps, attempt 
at eliminating them at each order through requiring final state clustering and looking for 
the optimal value of the jet-radius, and may be the jet veto too, such that non-global 
corrections are wiped out. In our rough approximation, we find that NGLs are completely 
eliminated, leaving only the primary Sudakov form factor, at R s < 3 x 10 -5 (equivalent 
to R < 0.01). Although this value is very small and not of any practical significance, 
including the full all-orders resummed results for both NGLs and CLs might result in 
practically larger values. Whether this is indeed the case remains to be investigated. If 
it turns out that the optimal radius is relatively large, 0.04 < R s (0.4 < R), then final 
state clustering will be the key to solve the NGLs subtlety of non-global observables. 

In the next section, we compare our analytical calculations to EVENT2. In particular, 
we focus on establishing the presence of NGLs and CLs in the te distribution at NLO. 

5 Numerical results 

The te numerical distribution has been computed using the fixed-order NLO QCD 
program EVENT2. The program implements the Catani-Seymour subtraction formalism 
for NLO corrections to two- and three-jet events observables in e + e~ annihilation. Final 
state partons have been clustered into jets using the Fast Jet library [39]. The latter 
provides an implementation of the longitudinally invariant k t , Cambridge-Aachen (CA) 
and anti— kt jet finders along with many others. Cone algorithms such as SISCone [40] are 
also implemented as plugins for the package. It should be noted that the e + e~ version 
of the aforementioned algorithms employs the following clustering condition for a pair of 
partons (ij) 

1 — cos% < 1 — cos(.R), (5.1) 

17 The discrepancy at R s — 0.30 is due to the fact that we have employed the small angles approximation 
in the C-A calculations. 
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Figure 5: Comparison of analytical resummed differential distribution d£/dr£ where: 
only primary term included (4.3), primary and NGLs factor included in the anti-A^ al- 
gorithm (4.3) and primary + NGLs + CLs factors included (4.6). The plots are shown 
for various jet-radii with a jet veto E = 0.1Q. The coupling is taken at the Z mass to 
be a s (Mz) — 0.118. The plots are only meant to give a rough estimate of the effects of 
NGLs in non-clustered as well as clustered final states. 



where R is the jet-radius parameter used in Fast Jet 18 . Compared to eqs (2.3) and (2.4), 
R = cos _1 (l — 2R S ). The exact numerical distributions (l/<7o)(d<j e /dL), with L = —L = 
ln(re ),. for the three colour channels, Cp,CpCA and CpT R Uf, have been obtained with 
10 11 events in the bin range > L > —14. We have used four values for the jet-radius: 
R s = 0.50, R s = 0.30, R s = 0.12 and R s = 0.04, with an energy veto E = 0.01 Q. 



Standard deviations on individual bins range from 10 4 % to 10 2 %. 

We plot the difference between the numerical and analytical distributions at both LO 
and NLO, 

r(L) = ^dT - ML' (5 ' 2) 

where da^/codL is given in eq. (B.6). Recall that at small values of the jet shape, r^ , 
the finite remainder function D^ n {TE ) is vanishingly small and will thus be ignored. For 
the case where the jet shape is global (R s = 0.50 and the threshold thrust reduces to 
thrust), we expect a full cancellation of singular terms and thus r should be a constant 
line corresponding to the NNLL coefficient (H21 in eq. (B.5)). For R s < 0.5, the jet shape 



18 In Fast Jet's manual R is allowed to go up to ir. Since we are interested in two-jet events the jet 
size cannot be wider than a hemisphere. Thus we restrict R to be less than ir/2. 
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Figure 6: The difference between EVENT2 and te LO distribution for various jet radii in 
both anti-fct (left) and CA (right) algorithms. 




Figure 7: The CpC& part of the difference between EVENT2 and (left) Te primary (global) 
distribution and (right) te distribution including NGLs for various jet radii in anti-fc* 
algorithm. 



is non-global and we expect r to have a slope if NGLs contribution is excluded. If our 
analytical calculations of the NGLs' coefficient, both for anti-fct and C-A algorithms, are 
correct then upon adding the latter to H 2 2 the slope should vanish and r becomes flat 
signalling a complete cancellation of terms up to NLL level. Similar behaviour should be 
seen with the CLs' coefficient C<f for the C-A algorithm case. Considering figs. 6 - 14, 
we make the following observations: 

• At LO, the distribution is independent of the jet definition. From eq. (B.5) we have 

H u = -3 + AL Rs = -4, -8.51, -14.75, -20.95; for R s = 0.5, 0.3, 0.12, 0.04. (5.3) 

Compared to the numerical results shown in fig. 6 we see a complete agreement. 
The cut-off in fig. 6 is due to the fact that at LO te < -R s /(1 + R s ) (eq. (2.10)). 

• For the NLO distribution in the anti-A^ algorithm, fig. 7 (left) illustrates the exis- 
tence of NGLs. The flatness of the r(L) curve, in fig. 7 (right), at L below about 
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R s Jet alg 


CpTnUf piece of 


CfCa piece of 


Cp piece of 




TT-analyt 


H 21 


TT-analyt 
n '2l 


rrnum 
H 2l 


Tj-analyt 
n 2l 




anti-fct 


5.20 ±0.14 


5.00 


-7.26 ±0.22 


-7.04 


-12.40 ±0.57 


-12.68 


C-A 


4.99 ±0.19 


-6.97 ±0.11 


-12.99 ±0.67 


0.30 


anti-fci 


13.92 ±0.15 


7.80 


-88.35 ±0.26 


-11.45 


62.73 ±0.51 


62.20 


C-A 


11.24 ±0.08 


-48.46 ±0.31 


76.38 ±0.50 


0.12 


anti-fcf 


15.75 ±0.12 


8.54 


-106.80 ±0.34 


-8.89 


384.04 ±0.90 


385.27 


C-A 


13.65 ±0.08 


-51.46 ±0.35 


405.48 ±0.57 


0.04 


anti-A^ 


12.86 ±0.19 


5.66 


-107.67 ±0.28 


3.59 


1017.93 ±0.53 


1022.43 


C-A 


10.60 ±0.23 


-45.67 ±0.23 


1044.33 ± 1.24 



Table 1: if 21 numerical vs analytical values for all three colour pieces. The numerical 
values were obtained through fitting the flat curve r(L) with the function (see eq. (B.7)) 
H™ m +e- l (B+C e~ l ). The analytical values, H^ alyt , are taken from eqs. (B.5) and (C.4) 
which only include the primary emission contribution with neither non-global nor clus- 
tering terms. R s = 0.50 corresponds to the global case and we expect the analytical and 
numerical values to be the same. 



—9 indicates a complete cancellation up to single log level. The Cp and CpT R rif 
pieces are shown in fig. 8. In table 1 we provide both numerical and analytical, 
taken from SCET calculations (C.4), values of the NNLL coefficient, H21, at the 
considered R s values for the three colour channels. It is evident from the table that 
there are sub leading i? s -dependent NGLs for both CpC^ and CpT R rif channels. 
Such logs have been analytically computed in [36] for the hemisphere jet mass. Our 
numerical results show that they are also present for finite-size jets 19 . The primary 
Cp channel is free from such subleading NGLs as numerical and analytical values 
of H21 in the anti-fct coincide. 

Notice that while the x-axis in all figures shown in this section corresponds to 
ln(T£ ) = log(rE ), i.e, the natural logarithm of the jet shape, that of [5] corre- 
sponds to the logarithm of base 10, log 10 (p) ~ log(p)/2, p = rp , of the jet shape. 
Given this, fig. 7 above is equivalent to fig. 7 of [5]. Neither Cp,CFT R rif plots nor 
subleading NGLs were considered in [5]. 

• The asymptotic region, i.e, the region where large logs are expected to dominate 
over non-logarithmic contributions, corresponds to L less than about —9 (for figs. 7, 
Cp piece in fig 8 and may even be less for the CpT R nf piece in fig. 8) and seems 
to decrease further as R s becomes smaller. A similar effect is seen in the thrust 
distribution, fig. 9, where the numerical distribution has been obtained using the 
full definition (2.5). 

• Considering the clustering case with C-A algorithm, figs. 10 and 11 illustrate the 
presence of CLs in the Cp channel. Clearly, the addition of CLs makes the remainder 
r flat in the region L < —9. To strengthen this observation even more, we plot in 

'as would be anticipated, since the finite-size jet mass is an extension to the hemisphere mass. 
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Figure 8: The (left) C F and (right) CpT R rif piece of the difference between EVENT2 and 
te distribution for various jet radii in the anti— kt algorithm. 



fig. 12 the difference between EVENT2 distributions in anti-A^ and C-A algorithms, 
^anti-kt /dL-dcr^ A /dL) /cr , for all colour pieces. The slopes for the Cf and C F C A 
indicate that an NLL positive i? s -dependent term, and possibly an NNLL term as 
well, have been induced by clustering. This is confirmed in table 1. Moreover, 
the fact that the difference between the latter distributions in the CfTrUj piece is 
non-vanishing implies an NNLL impact of clustering. 

Furthermore, we note from fig. 11 that C.f seems to slightly vary with the jet-radius 
parameter R s . This can be seen for large values of R s (R s = 0.3) where our small 
angles approximation (3.24) is not expected to apply. 

• Similar analysis to those carried in the anti— kt algorithm apply to the CpCa piece 
of the te distribution in the C-A algorithm. Including the NGLs makes the r(L) 
curve looks convincingly flat in the region L < —9, particularly for smaller values 
of R s , as shown in fig. 13. Recall that we have used the small R s limit in carrying 
out the computation of S%~ , eq. (3.32). Our findings agree with those reported 
in [5] for jet-radii up to R s ~ 0.3 (R ~ 1). 

For completeness, the CpT R nf piece of the r(L) in the C-A algorithm is depicted 
in fig. 14. As shown in table 1, clustering requirement again reduces the impact of 
the subleading NGLs in both CpCa arid CpTRrif channels. 

In summary, we have confirmed through explicit comparison to exact numerical distri- 
butions the existence of large NGLs and large CLs for the te distribution at NLL and 
beyond. In light of these findings, the surprising cancellation between primary-only an- 
alytical distribution and EVENT2 presented in [4] vl may be explained as follows. While 
each colour part [C F ,CpC a and CpTRrif) separately does not agree with EVENT2, as 
shown in figs. 10, 11, 13 and 14, their sum seems to agree with EVENT2 (recall that in 
of [4] vl only the sum of the three colour factors is plotted against EVENT2). Such an 
unexpected agreement can arise from the following possible sources: 
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Figure 9: The various colour pieces of the difference between EVENT2 and thrust distribu- 
tion using the full definition (2.5). The pQCD resummed analytical expression for thrust 
distribution can be found in, for example, [11]. 



• The ln(rE ) region considered in [4] vl does not correspond to the asymptotic region 
where large logs are expected to dominate over non-logarithmic terms. Thus the 
agreement shown in plot 1 of [4] vl does not convey any message and all one can 
say is that the non-logarithmic terms in the range [—9, 0] happen to cancel out (see 
fig 15). 

• NGLs are significantly reduced in the C — A algorithm especially for a jet radius 
R s = 0.3 (which is the one considered in [4] vl), as clearly seen in figs. 4 and 13 
(left). For smaller jet radii, there is a clear disagreement between the result of [4] 
vl and EVENT2 as shown in fig. 15 

We have also shown that clustering the final state partons with the C-A algorithm yielded 
a significant reduction in NGLs impact, at NLL and beyond, albeit inducing large CLs, 
at NLL and beyond, in the primary emission sector. 

6 Conclusion 

The jet mass with a jet veto, or simply the threshold thrust, is an example of a wider class 
of non-global observables. These have the characteristic of being sensitive to radiation 
into restricted regions of phase space, or sensitive to radiation into the whole phase space 



24 



(a s /2 7l) 2 C F 2 piece 



(a s /2 IE) C F piece 




Figure 10: The C F part of the difference between EVENT2 and (left) r^ primary (global) 
distribution and (right) r Eo distribution including CLs for various jet radii in the C-A 
algorithm. 
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Figure 11: Zoomed-in plots for the Cp part of the difference between EVENT2 and ana- 
lytical te distribution with and without CLs for (left) R s = 0.12 and (right) R s = 0.3 in 
the C-A algorithm. 

but differently in different regions. For such observables the universal Sudakov form factor 
fails to reproduce the full logarithmic structure even at NLL accuracy. New contributions 
that are dependent on various variables such as the jet size and jet definition appear at 
this logarithmic level. In this paper, we have elaborated on these very contributions for 
the aforementioned observable. 

Considering secondary emissions, we have computed the full analytical expression of 
the first term, S2, in a series of missing large logs, namely NGLs. The coefficient depends, 
as anticipated, on the jet size and saturates at its maximum in the limit where the 
latter, i.e, jet size, vanishes. This saturation value was used in [3] as an approximation 
to the full value in the small R s limit. It turns out that the approximation is valid 
for quite a wide range of R s . The formula for S2 has been checked against full exact 
numerical result obtained by the program EVENT2. The difference between the analytical 
and numerical differential distributions was shown to be asymptotically flat signalling a 
complete cancellation of singular terms up to NLL level. This has all been done for final 
states defined in the cone-like anti-fcj jet algorithm. 



25 




-8 -6 
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Figure 12: Plots of the three colour pieces of the difference between two EVENT2 distribu- 
tions corresponding to axiti-kt and C-A algorithms for various jet radii. We only show C F 
and CpT R rif results at three values of the jet-radius due to large errors in these colour 
channels. 

To illustrate the dependence of NLL on the jet definition, we have investigated the 
effects of applying the C-A algorithm on e + e~ final states. The impact of soft partons 
clustering is two-fold. On one side, it reduces the size of NGLs through shrinking the 
phase space region where the latter dominantly come from, i.e, the region where the 
emitter and emitted soft partons are just in and just out of the jet. On the other side, it 
gives rise to new NLL logarithmic contributions, CLs, in the primary emission sector. In 
the small jet-radius limit, the corresponding coefficient at second order has been shown, 
through comparison to EVENT2, to be independent of R s . 

Furthermore, our numerical analyses with EVENT2 have shown that the asymptotic 
region where the said large logs, in both anti-kt and C-A jet algorithms, dominate cor- 
responds to L < — 9 and decreases for smaller values of the jet-radius. As a by-product, 
we have found that there are subleading NGLs in both CfCa and CfTrUj pieces as well 
as subleading CLs in the C F piece of the te distribution. Clustering impact on NGLs 
has been observed to extend to NNLL level too. Regarding NGLs in CfCa channel in 
both jet algorithms, our findings serve as a confirmation of the corresponding calculations 
performed within SCET in [5]. 

Based on our rough approximation to NLL resummation, which is exponentiating the 
fixed-order result for both NGLs and CLs, it has been shown that it may be possible 
to completely eliminate the non-global correction to the primary Sudakov form factor 
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Figure 13: The CfCa part of the difference between EVENT2 and (left) t Eq primary (global) 
distribution and (right) r^ distribution including NGLs for various jet radii in the C-A 
algorithm. 
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Figure 14: The CfTrUj piece of the difference between EVENT2 and r^ distribution for 
various jet radii in C — A algorithm. 



at all-orders for events where final states clustering is applied. This elimination can 
be achieved by tuning the jet-radius parameter, R s , of the jet algorithm as well as the 
jet veto Eq. If such optimal values of R s and Eq are of practical significance, that is 
R ~ 0.4 or so and E ^> Aqcd, then the single-gluon exponentiation should be sufficient 
in describing the experimental data. A concrete answer of whether such optimal values 
exist can only be established once an all-orders resummation of primary, NGLs and CLs 
is performed. We postpone this investigation to future publications. 

As mentioned earlier in sec. 3 and shown in [3], the inclusive k t jet algorithm behaves 
in an identical way to C-A algorithm with regard to the threshold thrust distribution. It 
would be interesting to conduct similar studies for events defined in IRC cone algorithms 
such as the SISCone. In principle, one expects to see analogous effects not only for the 
threshold thrust but for all shape variables that are of non-global nature. Moreover, we 
reserve the extension of the findings of this paper to hadron-hadron collisions to future 
work. Apart from complications due to coloured initial state, we expect the gross features 
of this paper to apply. 
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A Derivation of LO distribution 



In the present section we outline the derivation of the full logarithmic part of the LO te 
integrated distribution (2.10). For the emission of a single gluon, i.e, e + e~ — > qqg, we 
define the kinematic variables, Xi = 2pi.Q/Q 2 = 2E i jQ and = 2pi.pj/Q 2 — 1 — x^ 
where k = l(q),2(q),3(g). The 0(a s ) matrix-element squared can be computed by 
considering two Feynman graphs corresponding to real emission of the gluon g off the 
two hard legs q, q. Applying the appropriate QCD Feynman rules and supplementing the 
three-body phase space factor, the corresponding differential distribution is given by 



dV 1 ) C F a s 



_[_ , y, 2 



crdxidxj 



2tt (l- Xl )(l-x 2 ) 



(a.i; 



where a is the total hadronic cross-section. Up to 0(a 2 ), it is given in terms of the Born 
cross-section, cr , by the relation [41] 



0o 



~2~ix 



3C f 
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with 



K 2 = -cf- + C F C A - 110 



K 2 + Q(al 



-y + 4Cs 



(A.2) 



(A.3) 
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The integration region, which is originally 1 > x±, x 2 > and x\ + x 2 > 1 and which leads 
to divergences, gets modified by introducing the jet shape variable. For three partons 
in the final state, te is zero unless two partons are clustered together. Therefore te 
is non-vanishing only in two-jet events. For the latter events, there are six ways of 
ordering the energy fractions Xi corresponding to six regions of phase space that needs to 
be integrated over. Due to X\ -H- x 2 symmetry of the matrix-element (A.l), one can only 
consider three regions and multiply the result by a factor of 2. These regions correspond 
to; Xi > x 2 > x 3 , X\ > £3 > x 2 and x 3 > X\ > x 2 . The threshold thrust is then given by 

te = (1 - £i)6(xi - x 2 , x 2 - £3)6 (2R S - 1 + cos 6 23 ) + 

+ (1 - xi)Q(xi - £ 3 , x 3 - x 2 )Q (2R S - 1 + cos #23) + 

+ (1 - £3)6(23 - Xl , x x - x 2 )Q {2R S - 1 + cos 12 ) , (A.4) 

where G(a— b, b—c) = Q(a—b)Q(b—c). To obtain the full logarithmic contribution it is suf- 
ficient to only consider regions where the gluon is the softest parton (x 3 = min(xj)). Other 
regions, last term in RHS of eq. (A.4), only contributes non-logarithmically. Adding up 
real and virtual contributions, in (2.8), one is only left with the virtual corrections in the 
range 0(1 — x\ — te )- The corresponding angular function in (A.4) may be written in 
terms of the energy fractions as, 

2(1 - Xl ) 2(1 - xi) c , 
1 - cos #23 = — « -i A.5 

x 2 x 3 x 3 

where the last approximation follows from the fact that the gluon is the softest, xi,x 2 ^> 
£3. Hence the two-jet contribution to the first order shape fraction is given by 

V m, P , C F a s f 1 -^ f Rs x l + x l n Rs 

Z w (t Eo , E ) = — / dx 2 / dxi — r6 — — - r E 

2vr Ji-r s (i-t Eo ) Ji+r Eo +x 2 (1 -x 1 ){l-x 2 ) \1 + R.s 

(A.6) 

In case of events with three-jets in the final state, the energy of the softest jet is 
vetoed to be less than Eq. The corresponding phase space constraint, left after real- 
virtual mis-cancellation, on the differential cross-section (A.l) reads 

Noting that x\ + x 2 + 23 = 2, one can obtain the corresponding integration limits on x\ 
and x 2 . Adding up the result of the latter integration with that of eq. (A.6) and making 
use of the following dilogarithm identities [42] 

Li 2 m + Li 2 (l — x) = — - hx(x) ln(l - x), 

o 

Li 2 (x) + Li 2 Q = ^-Iln 2 (x). (A.8) 

one obtains eq. (2.10). 
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B G nm coefficients 



The resultant coefficients from the expansion of the exponent in the resummed integrated 
distribution, eq. (4.5), are 
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(b.i; 



where L Rs = ln(i? s /(l — R s )) and = \il{2Eq/Q). The factor fo(R s ) only captures 
the first term of fo given in eq. (2.11). We simply replace fo H- fo when comparing to 
the numerical distribution. Moreover, we have introduced, for shorthand, the function 
So(R a ) given by (cf. eq. (3.12)), 



S 2 = —C f Ca S (R s ). 
The one-loop constant is given by, eq. (2.10), 



Ci 



c f 



■ 1 + y'' 



(B.2) 



(B.3) 



Expanding the total resummed distribution in eq. (4.5) to 0(a 2 ) and up to NLL we have 
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where (recall that L = ln(l/r^ ) => te = e L ) 

di(e-~ L ) 
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Differentiating (B.4) w.r.t. L, the NLO differential distribution reads 

%H = = «(Z) ft + (fi) A,(Z) + D fl (Z), 

dL cro dL v 2tt/ V2tt/ 

where the singular (logarithmic) terms are given by 

2 



Ds = 1 + (g) [Gl0+Cl] + (g 

^a(2) 

i> fl (Z) 



C20 + 7:^*10 + G1G10 + C2 



2Hi 2 L + Hii + ^di{e~ L ] 
dL 
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dL 



C Threshold thrust distribution in SCET 



(B.5) 



(B.6) 



(B.7) 



The resummation of the threshold thrust in SCET is presented in the current section 
for comparison with pQCD. We shall only present the final form of the resummed result 
taken from Refs. [4,24,27]. For a full derivation and more in depth discussion one should 
consult the latter references. The only task we have performed here is the expansion of 
the full resummed distribution to 0(a 2 s ). 
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C.l Resummation 

The general formula of the resummed distribution for the threshold thrust is given by [4, 
27] 



dZ SCET (T Eo ,R) _ da SCET 

OodT£ 



dr 



E 



exp [4S(fi h ,fij) + AS(ii s ,fij) - 4A H (n h ,[i s ) + AAj(fij, /i s )} 



R f 



-2 Ar (fluids) 



-2A r (p h ,Hj) 



x 



3 ln ^~ + d v^j 



§ r Eo (9 v ,fI s ) — 



TE Q\ v e-^ 



(C.l) 



See [4,27] for full notation. In order to compute the fixed-order expansion of (C.l) up 
to 0(a1), all scales should be set equal (/i^ = fij = fi s = Q). In this limit, the evolution 
factors S, Aj and Ah vanish. The differentiation w.r.t. 77 is carried out using the explicit 
form of j and . The final result of the integrated distribution may be cast in the 
generic form (4.5) with the constants and coefficients of the logs given by 



c 2 



c f 



+ C F T R n f 



" 2 \ 




+ t) 


1 


3jt 2 


7T 4 


8 + 




72 




7vr 


I 81 


6 



CfCa 



493 85vr 2 



324 



24 



1 



C in + c out 



2 5 



73?r 4 



28 



(C.2) 
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Considering primary emission, the only missing piece in the distribution is the two-loop 
constants in the soft function, namely C™ and . 



References 

[1] S. D. Ellis, A. Hornig, C. Lee, C. K. Vermilion, and J. R. Walsh Phys. Lett. B689 
(2010) 82-89, [arXiv: 0912. 0262]. 

[2] S. D. Ellis, C. K. Vermilion, J. R. Walsh, A. Hornig, and C. Lee JEEP 11 (2010) 
101, [arXiv: 1001.0014]. 

[3] A. Banfi, M. Dasgupta, K. Khelifa-Kerfa, and S. Marzani JEEP 08 (2010) 064, 
[arXiv: 1004.3483]. 

[4] R. Kelley, M. D. Schwartz, and H. X. Zhu arXiv: 1102.0561. 

[5] A. Hornig, C. Lee, J. R. Walsh, and S. Zuberi arXiv: 1110 . 0004. 

[6] H.-n. Li, Z. Li, and C.-R Yuan Phys. Rev. Lett. 107 (2011) 152001, 
[arXiv: 1107.4535]. 

[7] M. Dasgupta and G. P. Salam J. Phys. G30 (2004) R143, [hep-ph/0312283]. 
[8] M. Beneke Phys. Rept. 317 (1999) 1-142, [hep-ph/9807443]. 
[9] B. Webber hep-ph/9411384. 
[10] M. Beneke and V. M. Braun hep-ph/00 10208. 

[11] R. K. Ellis, W. J. Stirling, and B. R. Webber, QCD and Collider Physics. 
Cambridge monographs on particle physics, nuclear physics and cosmology. 
Cambridge University Press, 1996. 

[12] A. Abdesselam et al. Eur. Phys. J. C71 (2011) 1661, [arXiv : 1012 . 5412]. 

[13] M. Dasgupta and G. P. Salam Phys. Lett. B512 (2001) 323-330, 
[hep-ph/0104277]. 

[14] S. Catani, L. Trentadue, G. Turnock, and B. R. Webber Nucl. Phys. B407 (1993) 
3-42. 

[15] A. Banfi, G. P. Salam, and G. Zanderighi JEEP 01 (2002) 018, [hep-ph/0112156]. 



33 



[16] A. Banfi, G. P. Salam, and G. Zanderighi JEEP 1006 (2010) 038, 
[arXiv: 1001.4082]. 

[17] D. de Florian and M. Grazzini Nucl.Phys. B704 (2005) 387-403, 
[hep-ph/0407241]. 

[18] A. Gehrmann-De Ridder, T. Gehrmann, and E. W. N. Glover JEEP 09 (2005) 056, 
[hep-ph/0505111]. 

A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, and G. Heinrich JEEP 
12 (2007) 094, [arXiv:0711.471l]. 

[19] P. F. Monni, T. Gehrmann, and G. Luisoni JEEP 08 (2011) 010, 
[arXiv: 1105.4560]. 

[20] C. W. Bauer, S. Fleming, D. Pirjol, and I. W. Stewart Phys. Rev. D63 (2001) 
114020, [hep-ph/00 11336]. 

C. W. Bauer, D. Pirjol, and I. W. Stewart Phys. Rev. D65 (2002) 054022, 
[hep-ph/0109045]. 

[21] Y.-T. Chien and M. D. Schwartz JEEP 08 (2010) 058, [arXiv : 1005 . 1644]. 

R. Abbate, M. Fickinger, A. Hoang, V. Mateu, and I. W. Stewart PoS 
RADCOR2009 (2010) 040, [arXiv: 1004.4894]. 

[22] M. Cacciari and G. P. Salam Phys. Lett. B641 (2006) 57-61, [hep-ph/0512210]. 

[23] M. Cacciari, G. P. Salam, and G. Soyez JEEP 04 (2008) 063, [arXiv: 0802. 1189]. 

[24] R. Kelley and M. D. Schwartz Phys. Rev. D83 (2011) 033001, [arXiv: 1008.4355]. 

[25] S. D. Ellis and D. E. Soper Phys. Rev. D48 (1993) 3160-3166, [hep-ph/9305266]. 

[26] S. Catani and M. H. Seymour Phys. Lett. B378 (1996) 287-301, [hep-ph/9602277]. 

[27] T. Becher and M. D. Schwartz JEEP 07 (2008) 034, [arXiv : 0803 . 0342]. 

[28] A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, and G. Heinrich Phys. 
Rev. Lett. 99 (2007) 132002, [arXiv: 0707. 1285]. 

[29] G. Oderda and G. F. Sterman Phys. Rev. Lett. 81 (1998) 3591-3594, 
[hep-ph/9806530]. 

[30] M. Dasgupta and G. P. Salam JEEP 03 (2002) 017, [hep-ph/0203009]. 

[31] R. B. Appleby and M. H. Seymour JEEP 12 (2002) 063, [hep-ph/02 11426]. 

[32] A. Banfi and M. Dasgupta Phys. Lett. B628 (2005) 49-56, [hep-ph/0508159]. 

[33] J. C. Collins, D. E. Soper, and G. F. Sterman Nucl. Phys. B250 (1985) 199. 

[34] R. Bonciani, S. Catani, M. L. Mangano, and P. Nason Phys. Lett. B575 (2003) 
268-278, [hep-ph/0307035]. 

34 



[35] A. Banfi, G. Marchesini, and G. Smye JEEP 08 (2002) 006, [hep-ph/0206076]. 

[36] R. Kelley, M. D. Schwartz, R. M. Schabinger, and H. X. Zhu Phys. Rev. D84 
(2011) 045022, [arXiv: 1105.3676]. 

[37] A. Hornig, C. Lee, I. W. Stewart, J. R. Walsh, and S. Zuberi JEEP 08 (2011) 054, 
[arXiv: 1105.4628]. 

[38] Y. Delenda, R. Appleby, M. Dasgupta, and A. Banfi JEEP 0612 (2006) 044, 
[hep-ph/0610242]. 

[39] M. Cacciari, G. P. Salam, and G. Soyez arXiv : 1111 . 6097. 

[40] G. P. Salam and G. Soyez JEEP 05 (2007) 086, [arXiv: 0704. 0292]. 

[41] T. Appelquist and H. Georgi Phys. Rev. D8 (1973) 4000-4002. 

[42] http: //functions .wolfram. com/ZetaFunctionsandPolylogarithms/PolyLog2. 



35 



